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I.  STATEMENT  OF  THE  PROBLEM 

The  problem  here  under  investigation  can  be  simply  stated  as: 

"To  obtain  a  parametric  estimate  of  the  radar  cross  section  of  the 
ionized  wake  of  a  hypersonic  body  as  a  function  of  aspect  angle,  for  a  wide 
range  of  frequencies  and  values  of  the  physical  parameters  involved.  ” 

These  parameters  are: 

(a)  The  radar  frequency,  f  cps 

(b)  The  electron  collision  frequency,  v  sec  ^ 


(c)  The  electron  distribution  functions  and  their  parameters. 


Figure  1  Geometry  of  the  Problem 

It  will  be  assumed  that  the  wake  is  symmetric  about  the  z-axis  and 
therefore  at  a  small  distance  behind  the  body  (the  plane  z  =  o  in  Figure  1) 
the  electron  density  N  of  the  ionized  wake  can  be  described  functionedly  by 

N  =  N(z,  p,  t)  z  >  o 
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The  following  idealizations  will  be  now  imposed  on  the  problem: 

(i)  Quasi- static  approximation,  i.  e.  ,  take  N(  z,p,  t)  to  be  a  slowly- 
varying  function  of  time,  then 

N  ( z,  p,  t)  -*•  N  ( z,  p) 

(ii)  The  space  variation  of  N  is  assumed  to  be: 

N(z.  p)  =  N^f  (z)  ,^(p) 

and  for  preliminary  investigation,  let  f  and  0  take  the  form: 

z 

f  (z)  =  e'  ^  z  >  o  (1) 


and 


p  >  o 


(2) 


(iii)  The  ionization  may  be  replaced  by  the  macroscopic  complex 
dielectric  constant,  given  by: 


€  € 
o  r 


(p.  z)  =  e 


o 


N(p.  z)  ^ 
me  w(a)+  iv) 


cofto  +  iv )_ 


(3) 


wher'’ 


e 

o 


1 

36-11 


X  10 


€  (p.  z)  = 

r'r 

q  =  1.  6  X  10 
m  =  9- 1  X  10 


farad/meter,  is  the  permitivity  of 
free  space 

the  dielectric  constant  of  the  medium 

coulomb,  is  the  charge  per  electron 
kg,  is  the  mass  of  an  electron 
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CO 


rad/ sec 


IS  the  operating  radar  frequency 


(0 

P 


N(p. 


m  e 


rad/ sec 


is  the  plasma  frequency 


-i 

V  =  sec 


IS  the  collision  frequency 


With  the  above  assumptions,  the  parameters  involved  for  considera¬ 
tion  are: 


N  ,  a,  b,  CO  and  v 
o 

The  selection  of  Gaussian-Exponentail  dependence  for  the  electron  density 
obviously  is  not  the  only  possible  one  and  the  choice  was  based  simply  on 
the  fact  that  at  the  present  state  of  knowledge  it  seems  to  be  at  least  as 
realistic  as  any  other  and  that  it  yields  integrals  that  can  be  integrated  in 
closed  form.  However  more  complicated  parametric  representation, 
using  functions  f(z),  or  0  (p)  other  than  Equation  (2)  and  (3),  will  be  taken 
up  at  a  later  stage.  It  must  also  be  noted  that  due  to  the  variation  of  den¬ 
sity  with  altitude  one  should  include  an  altitude  dependence  of  the  electron 
density.  But  this  too  will  be  taken  into  consideration  later. 
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11.  MATHEMATICAL  FORMULATION 

Mathematically,  the  solution  to  the  problem  stated  in  (I)  is  the  solution 
to  Maxwell's  equations  for  an  inhomogeneous  medium.  Hence  a  brief  dis¬ 
cussion  of  the  problem  can  be  started  from  Maxwell's  Equations: 

V .  e  E  =  0  '] 


V<  H  =  0 


VxE  =  ioou  H 
—  f^o  — 


VxH=-'iwe€  E 

—  Q  r  — 


where  p  =  4-rr  x  1  0  is  the  permeability  of  free  space 


e  =  X  10  is  the  permittivity  of  free  space 


> 


II-  1 


e  =  The  dielectric  constant  of  the  medium,  and  the  assumptions  cr  =  0, 
r 

E  =  E  e  and  H  =  H  e  have  been  made. 


Elimination  of  E  (or  H)  in  the  above  equations  yields: 


VxVxE-k  E  =  0 


e  Vx  —  Vx  H  -  e  H  =  0 

r  €  —  r  — 

r 


II- 2 


4 


where  k 


2 

CO  u.  e  . 
o  o 


In  the  available  literature,  exact  solutions  for  the  above  set  of  equations 
are  given  only  for  simple  one -dimensional  problems.  However,  reasonably 
good  approximate  solutions  can  be  obtained  by  use  of  the  following  arguments: 


(a)  Few  algebraic  steps  reduce  II-2  into 


(1) 


Ve 


E  +  e  E  =  -  V  ■  E  ^ 

-  -/ 


?  ?  V  e 

V  H-fke  H= - - 

—  r  —  6 


X  Vx  H 


>  II-3 


Equation  II- 3  is  the  basis  of  geometric  optics  approximation.  The  argument 
is  that,  if 


Ve 

- -  «  k  II-4 

€ 

r 

(an  approximation  which  improves  as  cj  increases)  then  the  right  hand  side 
of  Equation  II-3  is  of  small  order  of  magnitude,  so  that  it  can  be  neglected. 

This  reduces  the  problem  into  a  solution  of  the  homogeneous  equettion: 

2  2 

V  0+ke  0~o 

where  the  standard  procedure  of  ray  tracing  may  apply. 

It  is  to  be  noted  that  in  most  applications  of  geometric  optics, 

>1,  therefore  Equation  II-4  can  be  satisfied  as  long  as  Ve^.  is  com¬ 
paratively  small.  In  the  equivalent  dielectric  constant  of  ionized  media, 
on  the  other  hand,  e  ^.  <1,  so  that  Equation  II-4  can  never  be  satisfied  near 
the  region  =  o,  even  when  Ve^  is  small.  In  so  far  as  we  know,  in  general, no 


Luneberg,  "Mathematical  Foundations  of  Geometric  optics,  "  Brown  Univer¬ 
sity  Lectures,  1944. 


5 


satisfactory  mathematical  discussion  exists  of  the  solution  of  Equation  II-3 

when  e  =  o,  but  the  physical  approximation  of  using  the  concept  of  critical 
r 

density  sometimes  yields  useful  numerical  solutions. 


(b)  An  alternate  approach,  more  amenable  to  numerical  calculations 
is  given  by  Khizhniak.  ^  This  is  done  by  rewriting  Equations  II- 2  into  the 
form: 


VxVxE-k^E  =k^(G^-l)  E 

VxVxH-k^H  =  -icoG  Vx{g  -1}  E 
—  —  o  r  — 


^  II- 5 


(From  here  on  for  notational  simplicity  g  implies  g^.) 

By  treating  the  right  hand  side  as  a  "source"  term,  the  standard  Green's 
function  technique  should  apply.  Using  the  free  space  Green's  function 


G{x,r') 


i  k(r  -  r') 
e  —  — 

4Tr  (r  r') 


II- 6 


where  r  and  r'  are  the  coordinates  of  the  source  and  point  of  observation 
respectively,.  Equations  II- 5  can  be  expressed  as  integral  equations: 


JE-  (£.)  =  ^ 


-4  [■&(£;)•  V'}V  G{_r, _£_')+  k^E(_r')G(_r,r_')]dV' 


II- 7 


H(_r)  =  H^(r)  -  i  (oG^  Vx  (£')  “  1]  E  iD 


Khizhniak,  N.  A.  ,  "Green's  Function  for  Maxwell's  Equations  for  In¬ 
homogeneous  Media.  "  Zhurnal  Tekhnicheskoi  Fiziki,7-28,  1958. 
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where  Eq  (r),  Hq  (r)  are  some  incident  field,  and  E  (r^),  H  (v)  are  the  total 
field.  If  we  introduce  the  notation  Eg(r)  and  (r^)  as  the  scattered  field, 
then  Equations  II-7  can  be  written  as  integral  equations  for  the  scattered 
field: 


-s‘-’  " 

V 

H  (r)  =  -iu)€ 


[{E^(r')+  E^  (£')}. V']V'G(r,£')+  k^{E^(r^)+ E^(r ')}G(r,  r -)  dV 

II-8 

'G{r,  r')x  {E^ir ')  +  E^(r ')}  [ e  U')-l]  dV 
V 


It  is  noted  that  equations  II-7  or  II- 8  have  a  very  clear  physical  interpreta¬ 
tion.  If  we  write, 

n  =  [e  (£')-!]  [E^(£;^)+  E^(£')]  G(£,£^)dV'  II-9 

V 


as  the  Hertzian  potential  due  to  induced  polarization 

then  it  is  obvious  that  Equations  II- 8  can  be  reduced  to 

E^  (r)  =  w-  n  +  n 

H  (r)  =-ico  €  Vx  n 
—  s  —  o  — 


II -10 


which  agrees  with  the  standard  formulation  of  the  radiation  problem. 

Equations  II- 8  are  integial  equations  of  the  Fredholm  type  for  the 
scattered  field  Eg  (r).  Perhaps  some  detailed  investigation  is  necessary 
for  the  nature  of  the  solutions  for  various  forms  of  [e  (r)  -Ij . 

For  immediate  numerical  answers,  we  found  that  this  formulation 
is  particularly  attractive  for  the  following  reasons: 
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(i)  For  the  region  where  je  (r)  -  1  j  <  1,  which  is  true  for  frequencies 
higher  than  the  critical  frequency,  approximate  solutions  can  be  obtained 
by  using  ^  (r)  in  place  of  E(r).  This  is  essentially  the  Born  approximation. 

(ii)  In  principle,  the  high  order  corrections  to  (i)  can  be  obtained  by 
standard  Neumann  series.  But  this  involves  too  much  numerical  work, 
unless  some  analytical  integration  can  be  done  first. 


(iii)  In  conjunction  with  (i),  which  actually  amounts  to  a  "guess"  of 
E(r)  in  the  integral  II-7  or  II-9,  a  reasonable  guess  of  E(r)  (or  more 
precisely  the  induced  polarization)  should  improve  the  results  greatly. 

A  classical  example  is  the  Kirchhoff  diffraction  theory  where  the  induced 
surface  current  on  a  body  is  guessed  from  the  incident  field.  A  reasonable 
"guess"  of  the  field  inside  the  ionized  region,  especially  at  the  region 
iwhere  6(r)  <  o,  will  be  proposed  later. 


(iv)  From  the  correct  guess  of  E,  the  scattering  cross  section  can 
be  obtained  by  evaluating  one  integral.  This  is  given  by; 


<r  (^'^)  = 


1 


X  CfC  [e(r')-l]  S{T^)e 


—  ik  ^  •  r  * 

£  (jV 


V 


II-ll 


where  §  '  is  the  direction  of  observation  and  E^  is  the  amplitude  of  the  in¬ 
cident  field. 


For  derivation  of  II-ll  see  Appendix. 
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III.  ESTIMATION  OF  THE  RADAR  CROSS-SECTION 

The  evaluation  of  the  radar  cross-section  (given  by  Equation  II- 11) 
depends  on  the  use  of  the  correct  expression  for  E(r).  As  a  first  order 
approximation  assume  that  for  z  >  o 

?.  <£)  =  Eo 

This  approximation  enables  one  to  evaluate  the  integral  in  Equation  (11-11) 
in  closed  form. 


Refer  now  to  Figure  2;  in  cylindrical  coordinates  (p,  <p,  z),  any  point 
in  space  is  represented  by 

r'  (p  cos  <p,  p  sin  <j>,  z) 

and  the  symbol  ~  should  be  read  "whose  components  are". 


The  space  z  >  o  is  then  replaced  by  a  pseudo  dielectric  for  which 

o  .2 


€(r>)  -  1  =  - 


N  q 
o 


me  a)(a)+iv) 
o 


a2 

e  e 


z 

b 


a2 

=  -  o:  e  ^  e 


z 

b 


where 


N  q 


o  ... 

a  =  — - — — — can  be  easily  calculated, 
me  w(wTiu) 
o 


Because  of  the  axial  symmetry  the  incident  plane  wave  may  be  taken  as 


(£') 


A 

E  e 


A 

e 


o 


ikr’  •  § 
e  — 

ik(p  cos  <j)  sin  0  +  z  cos  0  ) 

e  o  o 


9 


w> 


where  the  direction  of  the  incident  wave 


s  (sin  0  ,  0,  cos  0  ) 
o  o 

and  its  direction  of  polarization 

A 

e  cos  0  .  0,  sin  0  ) 
o  o 

are  assumed  to  be  in  the  x-z  plane. 

The  direction  of  observation  is  taken  as 


~  (cos  0  cos  <b,  sin  0  sin  <b  ,  cos  0  ) 
s  s  s  s  s 

Using  the  above  expressions,  the  volume  integral  in  Equation  II-ll, 
namely 


I  =  j'fy  [€  (r')  1]  E(r')  e 


A  »  » 

-IKS  ‘  r  ,  ? 

—  dv 


i«lf  spo.ce. 
Z  >  0 


can  be  easily  carried  out  in  steps.  Substituting 

_  _P ^ 

,  I.  ,  -2  b 

£  (r  )  ■  1  =  -  O'  e  ^  e 


E  (£)  =  E^  (£} 


_  A  i  [  k  z  cos  0  +  k  p  cos  6  sin  0  1 

it/  6  6  ^  O 

O 


and 


s  *  ‘  r'  =  p  sin  0  cos  ( <j)  ■  ‘j’  )  +  z  cos  0 
—  s  s  s 


2tt 


oo 


CO 


so  that, 

J[_  (s  ^  )  =  -E^  O'  e^d  <1)^  d  z  ^  p  d  p 


--  ^ 

a  b  ik  z  COS0O  ik  p  cos  4>  sin  0 
e  e  e  e  o 


o  o 


^  iii-i 


•  e 


-ikpsin0  cos  (  )  -ikzcos0 

s  s  s 


11 


1 


Integrating  over  z: 

CO  z 


I 


d  z  e  e 


ik  z  f  cos  6  -  cos  8  1  =  - 

^  o  s-'  1  . 


- ik  (cos  0  -cos  0 

b  os 


hence, 

I  A 


A  „ 

a  e  E 


2tt  « 


K^'.  s)  =  -- 


2  . 


/id  p  e  e 


ik  p  cos  (})  sin  0jj-sin  0gCOs((})-<|>g 


— -ik(cos0  -  cos  0  )o  o 
b  os 


To  integrate  over  <{>  one  writes: 

cos  <i)sin  0  -  sin  0  cos  ( cb-  <b  )  =  cos  Asin  0  -  sin  0  (  cos  icosi  +  sin  <t>sin  6  ) 

o  s  s  o  s  s  s 

=  cos  (b  f  sin  0  -  sin  0  cos  tb  1  -  sin  <b  f  6  sin  6  1 

=  d  COS  { tb  +  4*  J 

d 


where 


/  2  2  2 
d  =V  (sin  0  -  sin  0  cos  (b  )  t  sin  0  sin  <b 

o  s  s  s  s 


2  2 
sin  0-2  sin  0  sin  0  cos  <b  +  sin  0 
O  OSS  s 


and 


tan  (b 


sin  0  sin  4> 
s  s 


d  sin  0  -  sin  0  cos  <b 

o  s  s 


From  the  well  known  relation 
2it 

d  (b  e 


^  j  X  _i(kpd)  cos  ((b  +  «bj) 


d'  =  2  TT  J  (kdp) 
o 


one  obtains: 


00 


,  ,A ,  A  -a  e  2ir  Eo 

I  (s  ,  s  )  = 


I 


JL 


a 


—  -  ik  (cos  0  -cos  0 
b  o  so 


pdp  e  (kdp) 
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Finally,  integration  over  p  yields 


T  r  A  . 

I  (s  ,  s  )■ 


ae  2  IT  E 


o 


— -ik(cos0  -  cos  0  ) 
b  os 


III- 2 


Making  now  use  of  the  identity 


(s  *  X  e  )^  =  1  -  (s'*  *e  )^  =  1  -  [  cos  0  sin  0  -  cos  0  sin  0  cos  (j>  ] 


the  cross  section  can  be  written  as: 

2  -,Z 


A  ,  A 

0-  (s  ,  s  )  = 


Noq 


m  € 


2  2  2 
CO  (w  +v  ) 


/I  ^2  ,  2 ,  ^  -  ,2 

(—  )  +k  (cos  0  -  cos  0  ) 
b  os 


1  7 

__  e  ^ 
16 


k2d2a2 


III- 3 


{  1  -  [cos  0  sin  0  -  cos  0  sin  0  coscj)  ]  } 

^  c  O  OSS 


Equation  III-3  is  the  expression  we  were  looking  for.  It  is  the  first 
order  approximation  to  the  radar  cross  section  of  an  ionized  wake  viewed 
from  an  arbitrary  angle  0g  when  it  is  illuminated  by  an  electromagnetic 
plane  wave  incident  at  an  arbitrary  angle  0^. 

Higher  order  correction  terms  to  equation  III-3  will  be  the  object 
of  future  work. 


The  validity  of  Equation  III- 3  is  doubtful  for  all  ranges  of  parameters. 
It  is  true  if  the  following  conditions  hold: 

2  4-^ 


(a) 


m  e 


(b) 


CO  (cjj  +  i  V) 
V  <  <  CO 


<1 


III-4 


III-  5 


However,  on  account  of  its  simplicity  a  wide  range  of  calculations  should  be 
carried  out.  It  can  be  used  at  least  as  the  first  order  estimation  of  the 
aspect  effect  on  radar  cross-section. 
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In  the  c;ase  of  backs cattering,  which  is  of  special  interest, 


0  =  180°  -  0  and  d)  =  180° 

s  o  S 

Hence: 

cos  4>  =  -cos  0  ,  sin  0  =  sin  0  ,  cos  A  -  -1 

s  o  s  o  s 

d  =2  sin  0 

o 

tan  d) ,  =  0 
d 


and  assiiming  v  =  0,  the  cross  section  as  a  function  of  aspect  angle  is  given 
by: 


U) 


0-  (0  )  = 
o 


p 

16  c 


q_ 

4 


2  2  2 
- 2 a  k  sin  0 


— _  T  4k^  cos^  0 


III- 6 


where 


u)  r  rad/ sec,  the  plasma  frequency  at  z  =  0,  p  =  0 
^o 

8  , 

c  =  3  X  10  m/ sec,  the  velocity  of  light  in  vacuum 


c 


and 

(1) 

a  &  b  are  scale  parameters  of  the  electron  distribution. 

A  study  of  Equation  III-6  indicates  that  as  a  first  order  approxi¬ 
mation  the  following  conclusions  can  be  drawn: 


a.  The  cross  section  is  proportional  to  the  square  to  the  electron 
density  at  the  origin. 


(1) 


For  the  case  of  backs cattering  at  nose-on  incidence  (i.  e.  0  =0)  and  as¬ 
suming  b  >  >  1,  the  cross  section  reduces  to: 

1 

o-(O)  =—  III-7 


64 


"c'^ 


This  is  in  exact  agreement  with  the  result  obtained  by  S.  Altshuler  who  in¬ 
vestigated  this  special  case  starting  in  a  different  way. 


b.  The  cross  section  is  insensitive  to  the  wake  length,  except  at 
near  broadside  view  when  it  varies  as  the  square  of  the  length. 

c.  The  effect  of  frequency  on  the  cross  section  is  introduced  in 
two  ways,  by  an  inverse  square  relationship  and  by  an  ex¬ 
ponential  to  the  square  of  the  frequency.  Hence  the  cross 
section  decreases  rapidly  as  frequency  increases. 

d.  For  a  constant  flux  of  electrons  crossing  the  z  =  0  plane  (see 
next  section)  tOp  is  inversely  proportional  to  a.  Hence  the  ef¬ 
fect  of  a  enters  through  the  exponential  term  only,  and  the 
cross  sc  ':ion  decreases  as  a  increases. 


r 
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IV.  DETERMINATION  OF  BACKSCATTERING 

CROSS-SECTION 


For  demonstration  of  the  significance  of  the  obtained  results,  let 
us  substitute  values  for  the  parameters  in  equation  III- 6  and  calculate  the 
backscattering  cross  section  as  a  function  of  the  aspect  angle. 

There  are  four  parameters  that  must  be  specified  before  calculations 
can  be  made.  These  are:  f,  a,  and  b. 


To  obtain  a  reasonable  value  for  Nq  let  us  take  the  number  of 

7  1 

electrons  crossing  the  z  =  o  plane  to  be  10  ^  electrons/ sec,  a  value 
typical  of  wakes  of  hypersonic  bodies,  and  let  its  velocity  v  equal 
6,000  m/sec.  Then, 

■  "72 


V  N. 


=  10^^  el/ sec 


whicL.  after  integration  yields: 


No 


10 18 


el/m3 


IV- 1 


At  present,  flow  field  studies  of  wakes  indicate  that  reasonable 
values  for  a  and  b  are: 


0.2  <a  <1.  0  body  diameter 
5  ^b  -$'100  body  diameters 

and  an  even  greater  variation  is  the  value  of  b  ,  especially  at  the  high 
end,  is  possible. 

For  a  one-foot  diameter  sphere  these  ranges  translate  to: 

0.06  .$  a  $.0.30  meter 


1.5  $b  $30  meters 


Figures  3,  4,  5  and  6  show  the  variation  of  the  radar  cross  section 
as  a  functirii  of  aspect  angle  for  frequencies  in  the  UHF,  L,  C  and  X  band 
respectively,  with  a  as  parameter.  The  values  chosen  for  a  are: 

a  =  0.  06  m 
and  a  =  0.  15  m 

and  represent  the  range  of  most  probable  values  of  a. 

A  study  of  Figures  3  to  8  shows  that  the  cross  section  has  a  maxi¬ 
mum  at  Gq  =  0°,  then  it  decreases  as  Gq  increases,  reaches  a  minimum 
somewhere  between  45°  and  85°  (the  lower  the  frequency  the  earlier  the 
minimum  occurs)  then  increases  at  broadside  view.  However,  due  to  the 
exponential  relationship  to  the  square  of  the  frequency  the  cross  section 
falls  off  so  rapidly  with  aspect  angle  for  frequencies  higher  than  UHF  range 
that  it  is  practically  zero  for  Gq  10^  and  although  it  goes  up  at  Go  =  90° 
by  five  or  even  six  orders  of  magnitude,  depending  on  b,  the  cross  section 
remains  practically  zero  (e.g.  for  f  =  1.  5  KMC/S  and  a  =  0,  15  m,  £r(85°) 
m^  and  0'*(90°)  m^).  Hence  the  curves  in  Figures  3  to 

7  are  valid  for  any  of  the  possible  values  which  b  may  be  expected  to  take. 

That  the  backscattering  cross-section  should  have  a  relative  maxi¬ 
mum  for  Gq  =  O  was  expected  because  of  the  assumed  geometry  of  the 
wake,  having  a  plane  surface  with  the  highest  electron  density  at  the  plane 
z  =  O. 


However  what  was  unexpected  is  the  fact  that  the  cross-section  for 
0°  180°  is  symmetric  about  the  Gq  =  90°  value.  We  believe  that 

higher  order  corrections  on  the  cross-section  as  given  by  III-6  will  remove 
this  anomaly. 


F  or 

calculations 

into; 


Y  =■  o,  as  it  was  taken  to  be  in  the  backscattering  cross  section 
,  the  validity  conditions  given  by  Equat  ons  III- 4  and  III- 5  combine 


pZ 

Nq  e  ■ 

m  €o 

06  2 


e 


z 


c 


IV- 2 


If  this  inequality  is  to  be  satisfied  everywhere  in  the  half  space  z  ^  O, 
we  must  have: 


1 


IV- 3 


Nq  2 

m  £o  =  ^Po 


or 


IV-3a 


Inserting  in  equation  IV-3  the  value  calculated  for  Nq  from  equation  IV- 1 
and  taking  a  =  0.  15  m  one  obtains 


OJ 


Po 


=  8.66  X  10^®  rad/sec 


which  implies  that  for  these  values  of  Nq  a  one  must  use  in  equation 
III-5  frequencies  higher  than  13,8  KMC/S  for  the  resvilts  to  be  strictly 
valid.  However,  conditions  III-4  and  III-5  are  rather  stringent  and  the 
resxilts  are  fairly  good  for 

2 

‘^PO  2 

-ZF 


or  even  a  little  higher.  This  indicates  that  equation  III- 6  can  be  used  to 
obtain  radar  cross  sections  at  C  and  X  band  frequencies  with  fairly  good 
results . 

Currently  calculations  for  bistatic  cross-sections  are  being  carried  out 
and  will  be  reported  in  a  subsequent  Technical  Note. 
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APPENDIX 


Derivation  of  Equation  II- II. 


From  Equation  II- 9, 
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